Multi-omics bioactivity profile-based chemical grouping and read-across: a case study with Daphnia magna and azo dyes.

Grouping/read-across is widely used for predicting the toxicity of data-poor target substance(s) using data-rich source substance(s). While the chemical industry and the regulators recognise its benefits, registration dossiers are often rejected due to weak analogue/category justifications based largely on the structural similarity of source and target substances. Here we demonstrate how multi-omics measurements can improve confidence in grouping via a statistical assessment of the similarity of molecular effects. Six azo dyes provided a pool of potential source substances to predict long-term toxicity to aquatic invertebrates (Daphnia magna) for the dye Disperse Yellow 3 (DY3) as the target substance. First, we assessed the structural similarities of the dyes, generating a grouping hypothesis with DY3 and two Sudan dyes within one group. Daphnia magna were exposed acutely to equi-effective doses of all seven dyes (each at 3 doses and 3 time points), transcriptomics and metabolomics data were generated from 760 samples. Multi-omics bioactivity profile-based grouping uniquely revealed that Sudan 1 (S1) is the most suitable analogue for read-across to DY3. Mapping ToxPrint structural fingerprints of the dyes onto the bioactivity profile-based grouping indicated an aromatic alcohol moiety could be responsible for this bioactivity similarity. The long-term reproductive toxicity to aquatic invertebrates of DY3 was predicted from S1 (21-day NOEC, 40 µg/L). This prediction was confirmed experimentally by measuring the toxicity of DY3 in D. magna. While limitations of this 'omics approach are identified, the study illustrates an effective statistical approach for building chemical groups.

A novel method to derive a human safety limit for PFOA by gene expression profiling and modelling.

Perfluorooctanoic acid (PFOA) is a persistent environmental contaminant that can accumulate in the human body due to its long half-life. This substance has been associated with liver, pancreatic, testicular and breast cancers, liver steatosis and endocrine disruption. PFOA is a member of a large group of substances also known as "forever chemicals" and the vast majority of substances of this group lack toxicological data that would enable their effective risk assessment in terms of human health hazards. This study aimed to derive a health-based guidance value for PFOA intake (ng/kg BW/day) from in vitro transcriptomics data. To this end, we developed an in silico workflow comprising five components: (i) sourcing in vitro hepatic transcriptomics concentration-response data; (ii) deriving molecular points of departure using BMDExpress3 and performing pathway analysis using gene set enrichment analysis (GSEA) to identify the most sensitive molecular pathways to PFOA exposure; (iii) estimating freely-dissolved PFOA concentrations in vitro using a mass balance model; (iv) estimating in vivo doses by reverse dosimetry using a PBK model for PFOA as part of a quantitative in vitro to in vivo extrapolation (QIVIVE) algorithm; and (v) calculating a tolerable daily intake (TDI) for PFOA. Fourteen percent of interrogated genes exhibited in vitro concentration-response relationships. GSEA pathway enrichment analysis revealed that "fatty acid metabolism" was the most sensitive pathway to PFOA exposure. In vitro free PFOA concentrations were calculated to be 2.9% of the nominal applied concentrations, and these free concentrations were input into the QIVIVE workflow. Exposure doses for a virtual population of 3,000 individuals were estimated, from which a TDI of 0.15 ng/kg BW/day for PFOA was calculated using the benchmark dose modelling software, PROAST. This TDI is comparable to previously published values of 1.16, 0.69, and 0.86 ng/kg BW/day by the European Food Safety Authority. In conclusion, this study demonstrates the combined utility of an "omics"-derived molecular point of departure and in silico QIVIVE workflow for setting health-based guidance values in anticipation of the acceptance of in vitro concentration-response molecular measurements in chemical risk assessment.

Demonstrating the reliability of in vivo metabolomics based chemical grouping: towards best practice.

While grouping/read-across is widely used to fill data gaps, chemical registration dossiers are often rejected due to weak category justifications based on structural similarity only. Metabolomics provides a route to robust chemical categories via evidence of shared molecular effects across source and target substances. To gain international acceptance, this approach must demonstrate high reliability, and best-practice guidance is required. The MetAbolomics ring Trial for CHemical groupING (MATCHING), comprising six industrial, government and academic ring-trial partners, evaluated inter-laboratory reproducibility and worked towards best-practice. An independent team selected eight substances (WY-14643, 4-chloro-3-nitroaniline, 17α-methyl-testosterone, trenbolone, aniline, dichlorprop-p, 2-chloroaniline, fenofibrate); ring-trial partners were blinded to their identities and modes-of-action. Plasma samples were derived from 28-day rat tests (two doses per substance), aliquoted, and distributed to partners. Each partner applied their preferred liquid chromatography-mass spectrometry (LC-MS) metabolomics workflows to acquire, process, quality assess, statistically analyze and report their grouping results to the European Chemicals Agency, to ensure the blinding conditions of the ring trial. Five of six partners, whose metabolomics datasets passed quality control, correctly identified the grouping of eight test substances into three categories, for both male and female rats. Strikingly, this was achieved even though a range of metabolomics approaches were used. Through assessing intrastudy quality-control samples, the sixth partner observed high technical variation and was unable to group the substances. By comparing workflows, we conclude that some heterogeneity in metabolomics methods is not detrimental to consistent grouping, and that assessing data quality prior to grouping is essential. We recommend development of international guidance for quality-control acceptance criteria. This study demonstrates the reliability of metabolomics for chemical grouping and works towards best-practice.

Simultaneously discovering the fate and biochemical effects of pharmaceuticals through untargeted metabolomics.

Untargeted metabolomics is an established approach in toxicology for characterising endogenous metabolic responses to xenobiotic exposure. Detecting the xenobiotic and its biotransformation products as part of the metabolomics analysis provides an opportunity to simultaneously gain deep insights into its fate and metabolism, and to associate the internal relative dose directly with endogenous metabolic responses. This integration of untargeted exposure and response measurements into a single assay has yet to be fully demonstrated. Here we assemble a workflow to discover and analyse pharmaceutical-related measurements from routine untargeted UHPLC-MS metabolomics datasets, derived from in vivo (rat plasma and cardiac tissue, and human plasma) and in vitro (human cardiomyocytes) studies that were principally designed to investigate endogenous metabolic responses to drug exposure. Our findings clearly demonstrate how untargeted metabolomics can discover extensive biotransformation maps, temporally-changing relative systemic exposure, and direct associations of endogenous biochemical responses to the internal dose.

Applying genomics in regulatory toxicology: a report of the ECETOC workshop on omics threshold on non-adversity.

In a joint effort involving scientists from academia, industry and regulatory agencies, ECETOC's activities in Omics have led to conceptual proposals for: (1) A framework that assures data quality for reporting and inclusion of Omics data in regulatory assessments; and (2) an approach to robustly quantify these data, prior to interpretation for regulatory use. In continuation of these activities this workshop explored and identified areas of need to facilitate robust interpretation of such data in the context of deriving points of departure (POD) for risk assessment and determining an adverse change from normal variation. ECETOC was amongst the first to systematically explore the application of Omics methods, now incorporated into the group of methods known as New Approach Methodologies (NAMs), to regulatory toxicology. This support has been in the form of both projects (primarily with CEFIC/LRI) and workshops. Outputs have led to projects included in the workplan of the Extended Advisory Group on Molecular Screening and Toxicogenomics (EAGMST) group of the Organisation for Economic Co-operation and Development (OECD) and to the drafting of OECD Guidance Documents for Omics data reporting, with potentially more to follow on data transformation and interpretation. The current workshop was the last in a series of technical methods development workshops, with a sub-focus on the derivation of a POD from Omics data. Workshop presentations demonstrated that Omics data developed within robust frameworks for both scientific data generation and analysis can be used to derive a POD. The issue of noise in the data was discussed as an important consideration for identifying robust Omics changes and deriving a POD. Such variability or "noise" can comprise technical or biological variation within a dataset and should clearly be distinguished from homeostatic responses. Adverse outcome pathways (AOPs) were considered a useful framework on which to assemble Omics methods, and a number of case examples were presented in illustration of this point. What is apparent is that high dimension data will always be subject to varying processing pipelines and hence interpretation, depending on the context they are used in. Yet, they can provide valuable input for regulatory toxicology, with the pre-condition being robust methods for the collection and processing of data together with a comprehensive description how the data were interpreted, and conclusions reached.

A proposed framework to evaluate the quality and reliability of targeted metabolomics assays from the UK Consortium on Metabolic Phenotyping (MAP/UK).

Targeted metabolite assays that measure tens or hundreds of pre-selected metabolites, typically using liquid chromatography-mass spectrometry, are increasingly being developed and applied to metabolic phenotyping studies. These are used both as standalone phenotyping methods and for the validation of putative metabolic biomarkers obtained from untargeted metabolomics studies. However, there are no widely accepted standards in the scientific community for ensuring reliability of the development and validation of targeted metabolite assays (referred to here as 'targeted metabolomics'). Most current practices attempt to adopt, with modifications, the strict guidance provided by drug regulatory authorities for analytical methods designed largely for measuring drugs and other xenobiotic analytes. Here, the regulatory guidance provided by the European Medicines Agency, US Food and Drug Administration and International Council for Harmonisation of Technical Requirements for Pharmaceuticals for Human Use are summarized. In this Perspective, we have adapted these guidelines and propose a less onerous 'tiered' approach to evaluate the reliability of a wide range of metabolomics analyses, addressing the need for community-accepted, harmonized guidelines for tiers other than full validation. This 'fit-for-purpose' tiered approach comprises four levels-discovery, screening, qualification and validation-and is discussed in the context of a range of targeted and untargeted metabolomics assays. Issues arising with targeted multiplexed metabolomics assays, and how these might be addressed, are considered. Furthermore, guidance is provided to assist the community with selecting the appropriate degree of reliability for a series of well-defined applications of metabolomics.

Derivation of metabolic point of departure using high-throughput in vitro metabolomics: investigating the importance of sampling time points on benchmark concentration values in the HepaRG cell line.

Amongst omics technologies, metabolomics should have particular value in regulatory toxicology as the measurement of the molecular phenotype is the closest to traditional apical endpoints, whilst offering mechanistic insights into the biological perturbations. Despite this, the application of untargeted metabolomics for point-of-departure (POD) derivation via benchmark concentration (BMC) modelling is still a relatively unexplored area. In this study, a high-throughput workflow was applied to derive PODs associated with a chemical exposure by measuring the intracellular metabolome of the HepaRG cell line following treatment with one of four chemicals (aflatoxin B1, benzo[a]pyrene, cyclosporin A, or rotenone), each at seven concentrations (aflatoxin B1, benzo[a]pyrene, cyclosporin A: from 0.2048 µM to 50 µM; rotenone: from 0.04096 to 10 µM) and five sampling time points (2, 6, 12, 24 and 48 h). The study explored three approaches to derive PODs using benchmark concentration modelling applied to single features in the metabolomics datasets or annotated metabolites or lipids: (1) the 1st rank-ordered unannotated feature, (2) the 1st rank-ordered putatively annotated feature (using a recently developed HepaRG-specific library of polar metabolites and lipids), and (3) 25th rank-ordered feature, demonstrating that for three out of four chemical datasets all of these approaches led to relatively consistent BMC values, varying less than tenfold across the methods. In addition, using the 1st rank-ordered unannotated feature it was possible to investigate temporal trends in the datasets, which were shown to be chemical specific. Furthermore, a possible integration of metabolomics-driven POD derivation with the liver steatosis adverse outcome pathway (AOP) was demonstrated. The study highlights that advances in technologies enable application of in vitro metabolomics at scale; however, greater confidence in metabolite identification is required to ensure PODs are mechanistically anchored.

Data-driven learning of narcosis mode of action identifies a CNS transcriptional signature shared between whole organism Caenorhabditis elegans and a fish gill cell line.

With the large numbers of man-made chemicals produced and released in the environment, there is a need to provide assessments on their potential effects on environmental safety and human health. Current regulatory frameworks rely on a mix of both hazard and risk-based approaches to make safety decisions, but the large number of chemicals in commerce combined with an increased need to conduct assessments in the absence of animal testing makes this increasingly challenging. This challenge is catalysing the use of more mechanistic knowledge in safety assessment from both in silico and in vitro approaches in the hope that this will increase confidence in being able to identify modes of action (MoA) for the chemicals in question. Here we approach this challenge by testing whether a functional genomics approach in C. elegans and in a fish cell line can identify molecular mechanisms underlying the effects of narcotics, and the effects of more specific acting toxicants. We show that narcosis affects the expression of neuronal genes associated with CNS function in C. elegans and in a fish cell line. Overall, we believe that our study provides an important step in developing mechanistically relevant biomarkers which can be used to screen for hazards, and which prevent the need for repeated animal or cross-species comparisons for each new chemical.

Microbiome function predicts amphibian chytridiomycosis disease dynamics.

BACKGROUND: The fungal pathogen Batrachochytrium dendrobatidis (Bd) threatens amphibian biodiversity and ecosystem stability worldwide. Amphibian skin microbial community structure has been linked to the clinical outcome of Bd infections, yet its overall functional importance is poorly understood. METHODS: Microbiome taxonomic and functional profiles were assessed using high-throughput bacterial 16S rRNA and fungal ITS2 gene sequencing, bacterial shotgun metagenomics and skin mucosal metabolomics. We sampled 56 wild midwife toads (Alytes obstetricans) from montane populations exhibiting Bd epizootic or enzootic disease dynamics. In addition, to assess whether disease-specific microbiome profiles were linked to microbe-mediated protection or Bd-induced perturbation, we performed a laboratory Bd challenge experiment whereby 40 young adult A. obstetricans were exposed to Bd or a control sham infection. We measured temporal changes in the microbiome as well as functional profiles of Bd-exposed and control animals at peak infection. RESULTS: Microbiome community structure and function differed in wild populations based on infection history and in experimental control versus Bd-exposed animals. Bd exposure in the laboratory resulted in dynamic changes in microbiome community structure and functional differences, with infection clearance in all but one infected animal. Sphingobacterium, Stenotrophomonas and an unclassified Commamonadaceae were associated with wild epizootic dynamics and also had reduced abundance in laboratory Bd-exposed animals that cleared infection, indicating a negative association with Bd resistance. This was further supported by microbe-metabolite integration which identified functionally relevant taxa driving disease outcome, of which Sphingobacterium and Bd were most influential in wild epizootic dynamics. The strong correlation between microbial taxonomic community composition and skin metabolome in the laboratory and field is inconsistent with microbial functional redundancy, indicating that differences in microbial taxonomy drive functional variation. Shotgun metagenomic analyses support these findings, with similar disease-associated patterns in beta diversity. Analysis of differentially abundant bacterial genes and pathways indicated that bacterial environmental sensing and Bd resource competition are likely to be important in driving infection outcomes. CONCLUSIONS: Bd infection drives altered microbiome taxonomic and functional profiles across laboratory and field environments. Our application of multi-omics analyses in experimental and field settings robustly predicts Bd disease dynamics and identifies novel candidate biomarkers of infection. Video Abstract.

Spatially Mapping the Baseline and Bisphenol-A Exposed Daphnia magna Lipidome Using Desorption Electrospray Ionization-Mass Spectrometry.

Untargeted lipidomics has previously been applied to the study of daphnids and the discovery of biomarkers that are indicative of toxicity. Typically, liquid chromatography-mass spectrometry is used to measure the changes in lipid abundance in whole-body homogenates of daphnids, each only ca. 3 mm in length which limits any biochemical interpretation of site-specific toxicity. Here, we applied mass spectrometry imaging of Daphnia magna to combine untargeted lipidomics with spatial resolution to map the molecular perturbations to defined anatomical regions. A desorption electrospray ionization-mass spectrometry (DESI-MS) method was optimized and applied to tissue sections of daphnids exposed to bisphenol-A (BPA) compared to unexposed controls, generating an untargeted mass spectrum at each pixel (35 µm2/pixel) within each section. First, unique lipid profiles from distinct tissue types were identified in whole-body daphnids using principal component analysis, specifically distinguishing appendages, eggs, eye, and gut. Second, changes in the lipidome were mapped over four stages of normal egg development and then the effect of BPA exposure on the egg lipidome was characterized. The primary perturbations to the lipidome were annotated as triacylglycerides and phosphatidylcholine, and the distributions of the individual lipid species within these classes were visualized in whole-body D. magna sections as ion images. Using an optimized DESI-MS workflow, the first ion images of D. magna tissue sections were generated, mapping both their baseline and BPA-perturbed lipidomes.

Automated Sample Preparation and Data Collection Workflow for High-Throughput In Vitro Metabolomics.

Regulatory bodies have started to recognise the value of in vitro screening and metabolomics as two types of new approach methodologies (NAMs) for chemical risk assessments, yet few high-throughput in vitro toxicometabolomics studies have been reported. A significant challenge is to implement automated sample preparation of the low biomass samples typically used for in vitro screening. Building on previous work, we have developed, characterised and demonstrated an automated sample preparation and analysis workflow for in vitro metabolomics of HepaRG cells in 96-well microplates using a Biomek i7 Hybrid Workstation (Beckman Coulter) and Orbitrap Elite (Thermo Scientific) high-resolution nanoelectrospray direct infusion mass spectrometry (nESI-DIMS), across polar metabolites and lipids. The experimental conditions evaluated included the day of metabolite extraction, order of extraction of samples in 96-well microplates, position of the 96-well microplate on the instrument's deck and well location within a microplate. By using the median relative standard deviation (mRSD (%)) of spectral features, we have demonstrated good repeatability of the workflow (final mRSD < 30%) with a low percentage of features outside the threshold applied for statistical analysis. To improve the quality of the automated workflow further, small method modifications were made and then applied to a large cohort study (4860 sample infusions across three nESI-DIMS assays), which confirmed very high repeatability of the whole workflow from cell culturing to metabolite measurements, whilst providing a significant improvement in sample throughput. It is envisioned that the automated in vitro metabolomics workflow will help to advance the application of metabolomics (as a part of NAMs) in chemical safety, primarily as an approach for high throughput screening and prioritisation.

Knowledge-Driven Approaches to Create the MTox700+ Metabolite Panel for Predicting Toxicity.

Endogenous metabolite levels describe the molecular phenotype that is most downstream from chemical exposure. Consequently, quantitative changes in metabolite levels have the potential to predict mode-of-action and adversity, with regulatory toxicology predicated on the latter. However, toxicity-related metabolic biomarker resources remain highly fragmented and incomplete. Although development of the S1500+ gene biomarker panel has accelerated the application of transcriptomics to toxicology, a similar initiative for metabolic biomarkers is lacking. Our aim was to define a publicly available metabolic biomarker panel, equivalent to S1500+, capable of predicting pathway perturbations and/or adverse outcomes. We conducted a systematic review of multiple toxicological resources, yielding 189 proposed metabolic biomarkers from existing assays (BASF, Bowes-44, and Tox21), 342 biomarkers from databases (Adverse Outcome Pathway Wiki, Comparative Toxicogenomics Database, QIAGEN Ingenuity Pathway Analysis, and Toxin and Toxin-Target Database), and 435 biomarkers from the literature. Evidence mapping across all 8 resources generated a panel of 722 metabolic biomarkers for toxicology (MTox700+), of which 462 (64%) are associated with molecular pathways and 575 (80%) with adverse outcomes. Comparing MTox700+ and S1500+ revealed that 418 (58%) metabolic biomarkers associate with pathways shared across both panels, with further metabolites mapping to unique pathways. Metabolite reference standards are commercially available for 646 (90%) of the panel metabolites, and assays exist for 578 (80%) of these biomarkers. This study has generated a publicly available metabolic biomarker panel for toxicology, which through its future laboratory deployment, is intended to help build foundational knowledge to support the generation of molecular mechanistic data for chemical hazard assessment.

Integrating in vitro metabolomics with a 96-well high-throughput screening platform.

INTRODUCTION: High-throughput screening (HTS) is emerging as an approach to support decision-making in chemical safety assessments. In parallel, in vitro metabolomics is a promising approach that can help accelerate the transition from animal models to high-throughput cell-based models in toxicity testing. OBJECTIVE: In this study we establish and evaluate a high-throughput metabolomics workflow that is compatible with a 96-well HTS platform employing 50,000 hepatocytes of HepaRG per well. METHODS: Low biomass cell samples were extracted for metabolomics analyses using a newly established semi-automated protocol, and the intracellular metabolites were analysed using a high-resolution spectral-stitching nanoelectrospray direct infusion mass spectrometry (nESI-DIMS) method that was modified for low sample biomass. RESULTS: The method was assessed with respect to sensitivity and repeatability of the entire workflow from cell culturing and sampling to measurement of the metabolic phenotype, demonstrating sufficient sensitivity (> 3000 features in hepatocyte extracts) and intra- and inter-plate repeatability for polar nESI-DIMS assays (median relative standard deviation < 30%). The assays were employed for a proof-of-principle toxicological study with a model toxicant, cadmium chloride, revealing changes in the metabolome across five sampling times in the 48-h exposure period. To allow the option for lipidomics analyses, the solvent system was extended by establishing separate extraction methods for polar metabolites and lipids. CONCLUSIONS: Experimental, analytical and informatics workflows reported here met pre-defined criteria in terms of sensitivity, repeatability and ability to detect metabolome changes induced by a toxicant and are ready for application in metabolomics-driven toxicity testing to complement HTS assays.

An integrated host-microbiome response to atrazine exposure mediates toxicity in Drosophila.

The gut microbiome produces vitamins, nutrients, and neurotransmitters, and helps to modulate the host immune system-and also plays a major role in the metabolism of many exogenous compounds, including drugs and chemical toxicants. However, the extent to which specific microbial species or communities modulate hazard upon exposure to chemicals remains largely opaque. Focusing on the effects of collateral dietary exposure to the widely used herbicide atrazine, we applied integrated omics and phenotypic screening to assess the role of the gut microbiome in modulating host resilience in Drosophila melanogaster. Transcriptional and metabolic responses to these compounds are sex-specific and depend strongly on the presence of the commensal microbiome. Sequencing the genomes of all abundant microbes in the fly gut revealed an enzymatic pathway responsible for atrazine detoxification unique to Acetobacter tropicalis. We find that Acetobacter tropicalis alone, in gnotobiotic animals, is sufficient to rescue increased atrazine toxicity to wild-type, conventionally reared levels. This work points toward the derivation of biotic strategies to improve host resilience to environmental chemical exposures, and illustrates the power of integrative omics to identify pathways responsible for adverse health outcomes.

Energy Starvation in Daphnia magna from Exposure to a Lithium Cobalt Oxide Nanomaterial.

Growing evidence across organisms points to altered energy metabolism as an adverse outcome of metal oxide nanomaterial toxicity, with a mechanism of toxicity potentially related to the redox chemistry of processes involved in energy production. Despite this evidence, the significance of this mechanism has gone unrecognized in nanotoxicology due to the field's focus on oxidative stress as a universal─but nonspecific─nanotoxicity mechanism. To further explore metabolic impacts, we determined lithium cobalt oxide's (LCO's) effects on these pathways in the model organism Daphnia magna through global gene-expression analysis using RNA-Seq and untargeted metabolomics by direct-injection mass spectrometry. Our results show that a sublethal 1 mg/L 48 h exposure of D. magna to LCO nanosheets causes significant impacts on metabolic pathways versus untreated controls, while exposure to ions released over 48 h does not. Specifically, transcriptomic analysis using DAVID indicated significant enrichment (Benjamini-adjusted p ≤0.0.5) in LCO-exposed animals for changes in pathways involved in the cellular response to starvation (25 genes), mitochondrial function (70 genes), ATP-binding (70 genes), oxidative phosphorylation (53 genes), NADH dehydrogenase activity (12 genes), and protein biosynthesis (40 genes). Metabolomic analysis using MetaboAnalyst indicated significant enrichment (γ-adjusted p <0.1) for changes in amino acid metabolism (19 metabolites) and starch, sucrose, and galactose metabolism (7 metabolites). Overlap of significantly impacted pathways by RNA-Seq and metabolomics suggests amino acid breakdown and increased sugar import for energy production. Results indicate that LCO-exposed Daphnia respond to energy starvation by altering metabolic pathways, both at the gene expression and metabolite levels. These results support altered energy production as a sensitive nanotoxicity adverse outcome for LCO exposure and suggest negative impacts on energy metabolism as an important avenue for future studies of nanotoxicity, including for other biological systems and for metal oxide nanomaterials more broadly.

An Extensive Metabolomics Workflow to Discover Cardiotoxin-Induced Molecular Perturbations in Microtissues.

Discovering modes of action and predictive biomarkers of drug-induced structural cardiotoxicity offers the potential to improve cardiac safety assessment of lead compounds and enhance preclinical to clinical translation during drug development. Cardiac microtissues are a promising, physiologically relevant, in vitro model, each composed of ca. 500 cells. While untargeted metabolomics is capable of generating hypotheses on toxicological modes of action and discovering metabolic biomarkers, applying this technology to low-biomass microtissues in suspension is experimentally challenging. Thus, we first evaluated a filtration-based approach for harvesting microtissues and assessed the sensitivity and reproducibility of nanoelectrospray direct infusion mass spectrometry (nESI-DIMS) measurements of intracellular extracts, revealing samples consisting of 28 pooled microtissues, harvested by filtration, are suitable for profiling the intracellular metabolome and lipidome. Subsequently, an extensive workflow combining nESI-DIMS untargeted metabolomics and lipidomics of intracellular extracts with ultra-high performance liquid chromatography-mass spectrometry (UHPLC-MS/MS) analysis of spent culture medium, to profile the metabolic footprint and quantify drug exposure concentrations, was implemented. Using the synthetic drug and model cardiotoxin sunitinib, time-resolved metabolic and lipid perturbations in cardiac microtissues were investigated, providing valuable data for generating hypotheses on toxicological modes of action and identifying putative biomarkers such as disruption of purine metabolism and perturbation of polyunsaturated fatty acid levels.

Progress towards an OECD reporting framework for transcriptomics and metabolomics in regulatory toxicology.

Omics methodologies are widely used in toxicological research to understand modes and mechanisms of toxicity. Increasingly, these methodologies are being applied to questions of regulatory interest such as molecular point-of-departure derivation and chemical grouping/read-across. Despite its value, widespread regulatory acceptance of omics data has not yet occurred. Barriers to the routine application of omics data in regulatory decision making have been: 1) lack of transparency for data processing methods used to convert raw data into an interpretable list of observations; and 2) lack of standardization in reporting to ensure that omics data, associated metadata and the methodologies used to generate results are available for review by stakeholders, including regulators. Thus, in 2017, the Organisation for Economic Co-operation and Development (OECD) Extended Advisory Group on Molecular Screening and Toxicogenomics (EAGMST) launched a project to develop guidance for the reporting of omics data aimed at fostering further regulatory use. Here, we report on the ongoing development of the first formal reporting framework describing the processing and analysis of both transcriptomic and metabolomic data for regulatory toxicology. We introduce the modular structure, content, harmonization and strategy for trialling this reporting framework prior to its publication by the OECD.

Modeling the metabolic profile of Mytilus edulis reveals molecular signatures linked to gonadal development, sex and environmental site.

The monitoring of anthropogenic chemicals in the aquatic environment including their potential effects on aquatic organisms, is important for protecting life under water, a key sustainable development goal. In parallel with monitoring the concentrations of chemicals of concern, sentinel species are often used to investigate the biological effects of contaminants. Among these, bivalve molluscs such as mussels are filter-feeding and sessile, hence an excellent model system for measuring localized pollution. This study investigates the relationship between the metabolic state of the blue mussel (Mytilus edulis) and its physiology in different environments. We developed a computational model based on a reference site (relatively unpolluted) and integrated seasonal dynamics of metabolite relative concentrations with key physiological indicators and environmental parameters. The analysis of the model revealed that changes in metabolite levels during an annual cycle are influenced by water temperature and are linked to gonadal development. This work supports the importance of data-driven biology and its potential in environmental monitoring.

Acoustic Mist Ionization Mass Spectrometry for Ultrahigh-Throughput Metabolomics Screening.

Incorporating safety data early in the drug discovery pipeline is key to reducing costly lead candidate failures. For a single drug development project, we estimate that several thousand samples per day require screening (<10 s per acquisition). While chromatography-based metabolomics has proven value at predicting toxicity from metabolic biomarker profiles, it lacks sufficiently high sample throughput. Acoustic mist ionization mass spectrometry (AMI-MS) is an atmospheric pressure ionization approach that can measure metabolites directly from 384-well plates with unparalleled speed. We sought to implement a signal processing and data analysis workflow to produce high-quality AMI-MS metabolomics data and to demonstrate its application to drug safety screening. An existing direct infusion mass spectrometry workflow was adapted, extended, optimized, and tested, utilizing three AMI-MS data sets acquired from technical and biological replicates of metabolite standards and HepG2 cell lysates and a toxicity study. Driven by criteria to minimize variance and maximize feature counts, an algorithm to extract the pulsed scan data was designed; parameters for signal-to-noise-ratio, replicate filter, sample filter, missing value filter, and RSD filter were all optimized; normalization and batch correction strategies were adapted; and cell phenotype filtering was implemented to exclude high cytotoxicity samples. The workflow was demonstrated using a highly replicated HepG2 toxicity data set, comprising 2772 samples from exposures to 16 drugs across 9 concentrations and generated in under 5 h, revealing metabolic phenotypes and individual metabolite changes that characterize specific modes of action. This AMI-MS workflow opens the door to ultrahigh-throughput metabolomics screening, increasing the rate of sample analysis by approximately 2 orders of magnitude.

Environmentally Relevant Iron Oxide Nanoparticles Produce Limited Acute Pulmonary Effects in Rats at Realistic Exposure Levels.

Iron is typically the dominant metal in the ultrafine fraction of airborne particulate matter. Various studies have investigated the toxicity of inhaled nano-sized iron oxide particles (FeOxNPs) but their results have been contradictory, with some indicating no or minor effects and others finding effects including oxidative stress and inflammation. Most studies, however, did not use materials reflecting the characteristics of FeOxNPs present in the environment. We, therefore, analysed the potential toxicity of FeOxNPs of different forms (Fe3O4, α-Fe2O3 and γ-Fe2O3) reflecting the characteristics of high iron content nano-sized particles sampled from the environment, both individually and in a mixture (FeOx-mix). A preliminary in vitro study indicated Fe3O4 and FeOx-mix were more cytotoxic than either form of Fe2O3 in human bronchial epithelial cells (BEAS-2B). Follow-up in vitro (0.003, 0.03, 0.3 µg/mL, 24 h) and in vivo (Sprague-Dawley rats, nose-only exposure, 50 µg/m3 and 500 µg/m3, 3 h/d × 3 d) studies therefore focused on these materials. Experiments in vitro explored responses at the molecular level via multi-omics analyses at concentrations below those at which significant cytotoxicity was evident to avoid detection of responses secondary to toxicity. Inhalation experiments used aerosol concentrations chosen to produce similar levels of particle deposition on the airway surface as were delivered in vitro. These were markedly higher than environmental concentrations. No clinical signs of toxicity were seen nor effects on BALF cell counts or LDH levels. There were also no significant changes in transcriptomic or metabolomic responses in lung or BEAS-2B cells to suggest adverse effects.